Optimal estimation of losses at the ultimate quantum limit with non-Gaussian states 
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We address the estimation of the loss parameter of a bosonic channel probed by arbitrary signals. Unlike the 
optimal Gaussian probes, which can attain the ultimate bound on precision asymptotically either for very small 
or very large losses, we prove that Fock states at any fixed photon number saturate the bound unconditionally 
for any value of the loss. In the relevant regime of low-energy probes, we demonstrate that superpositions of the 
first low-lying Fock states yield an absolute improvement over any Gaussian probe. Such few-photon states can 
be recast quite generally as truncations of de-Gaussified photon subtracted states. 



PACS numbers: 03.65.Ta, 03.67.Hk, 42.50.Dv 

Introduction. — Suppose an experimenter is given different 
stations connected by channels manufactured before he/she 
joined the laboratory. Assume that he/she is allowed to 
routinely transmit light beams prepared e.g. in coherent, 
squeezed, or number states through these channels with 
the purpose of implementing some communication network. 
Given the state at one station, one finds out that in general 
the state has been altered at the next node. The problem is 
then to determine what sort of noise is affecting the transmis- 
sions. This is a typical issue of quantum parameter estima- 
tion, whose solution is clearly of direct interest to practical 
situations like the one described. The experimenter wishes to 
determine the optimal probe state that has to be sent through 
the channel, and the optimal measurement that needs to be 
performed at the output in order to estimate (after repeating 
the process N times) the loss parameter with the maximum 
possible precision. In the case of amplitude damping bosonic 
channels, Monras and Paris provided a solution to the problem 
in the particular case of Gaussian input probe states (displaced 
and squeezed vacuum states) 1 1]. Gaussian states are easier to 
engineer by quantum optical means than more sophisticated 
non-Gaussian states. On the other hand, the optimal mea- 
surement needed for loss estimation involves manipulations 
such as displacement and squeezing of the output signal, and 
photon-counting which is a non-Gaussian measurement not 
belonging to the standard toolbox of linear optics. At a fun- 
damental level, using Gaussian probe inputs allows to satu- 
rate the ultimate bound on precision only asymptotically in 
the unphysical limits of infinitesimal or infinite losses, while 
in the realistic regime of intermediate loss the Gaussian-based 
estimation is clearly suboptimal |1]. Therefore considering 
Gaussian inputs does not solve the important problem of opti- 
mal estimation of loss in bosonic channels, both on theoretical 
and practical grounds. 

In this work we study the estimation of loss in bosonic 
channels probed by arbitrary non-Gaussian states. For any 
energy of the probes, we show that there exist non-Gaussian 



states improving the estimation compared to Gaussian states 
in all regimes of loss. Specifically, we prove that Fock states 
\n) (which can be produced deterministic ally in the labora- 
tory |2]) are the truly optimal probes that attain the ultimate 
quantum limit exactly, for any n and any value of the loss. 
The optimal estimation then requires only photon-counting, 
resulting in a technological simplification compared to the 
Gaussian case. For low-energy probes (mean photon num- 
ber smaller than 1) we construct optimal superpositions of 
the first k low-lying Fock states which improve the estima- 
tion over the Gaussian case already for fc = 2, 3 and approach 
the ultimate limit in a much broader range of losses. Inter- 
estingly, we find that the optimal superpositions for k = 2 
coiTespond to qutrit-like, two-photon truncations of photon- 
subtracted states, showing that de-Gaussification procedures 
generally allow enhanced performance in the task of loss es- 
timation in quantum channels. This result adds to the di- 
verse existing instances of non-Gaussianity as a "power-up" 
for quantum information, encountered in the optimal cloning 
of coherent states yD, continuous variable teleportation j^, 
nonlocality tests ISO, and entanglement distillation l^. 

Bosonic channels and quantum estimation. — We consider 
a bosonic channel described by the master equation dp/dt = 
(7/2)£[a]p for quantum states p, where C[a]p — 2apa^ — 
a^ap — pa^a, a being the annihilation operator on the Fock 
space of a single bosonic mode. The aim of our study is the 
optimal estimation of the loss parameter 7, or equivalently of 
4> € (0,7r/2) defined by tan^ (f> = exp(7t) — 1. In terms 
of (f), the master equation reads dp/dcfi = tan (f>C[a]p, whose 
general solution is of the form 






-(cos0)'^^'^a"po(a^)"(cos0)'^^''. (1) 



Let us recall the basic elements of quantum estimation the- 
ory OTL ISp, and the relevant tools of interest for the present 
case llllsl]. The optimal estimation of is achieved asymptot- 
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FIG. 1: (color online) Quantum Fisher information ii((^) versus the loss parameter (/> for different values of the input energy n. Solid black 
line: optimal Gaussian probes [ 1]; Thin gray line: coherent states; Dashed (blue) line: Fock states |n); Dotted (magenta) line: qubit-like states 
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I'i/'q ); Dot-dashed (orange) line: optimal qutrit-like states \ipQ ); Dot-dot-dashed (cyan) line: optimal quartet-like states |i/)q ) 
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ically by sending N independent and identically distributed 
optimal probe states po into the channel, and performing at 
each run the optimal measurement on the output signals p^, 
in order to construct an estimator to infer the true value of 
(/) with minimal variance. For any given po, the optimal out- 
put measurement can be exactly determined in terms of the 
symmetric logarithmic derivative (SLD) A((/)), defined im- 
plicitly as the Hermitian operator that satisfies dp^/d^ = 
(l/2)[p0A(0) + A(0)/90]. Using the spectral decomposition 
P4> = I]fc Pk\'4'k){i'k I one finds for the SLD: 



A((/)) = 2tan(/)^ 



pq 



{il)q\C[a]p\^p) 

Pp + Pq 



IV'9)(V'p| 



(2) 



The resulting minimum variance saturates the quantum 
Cramer-Rao bound, Var0[0] > l/[NH{4i)], where the 
quantum Fisher information (QFI) -/?((/)) reads H{4)) — 
Tr[(O0A((/))^]. The problem is thus recast in the determination 
of the optimal single-mode pure input states with a given finite 
mean energy (or mean photon number n), such that the QFI 
of the corresponding output states is maximal. The ultimate 
quantum limit on precision, that is computable in the ideal as- 
sumption that the experimenter may have access also to the 
degrees of freedom of the environment (i.e. to the oscillators 
internal to the channel), is achieved for estimators with [ 1 ] 
Var0[0] > l/(4nA^). This means that a truly optimal estima- 
tion requires input probes which yield at the output a QFI ex- 
actly equal to 4n (for any single run). If the ensemble of input 
signals is limited to Gaussian states 1 1 ], the estimation is never 
optimal: the ultimate limit is attained only asymptotically for 
(^ approaching or 7r/2, while i?((/)) for the best Gaussian 
probes can get as low as ^ 2n for intermediate losses (see 
Fig.[T]i. Here we show that non-Gaussian probes, Fock states 
and low-lying superpositions thereof, are indeed optimal for 
the estimation of loss in bosonic channels with the maximum 
precision allowed by laws of quantum mechanics. 

Fock states. — Let us consider, as input probes, 
Fock states po = l"-)("-l (^^ = ''^)- The 



evolved state, according to Eq. ([T), reads p^ = 
ELo(sin^'?^)''(fc)(cos2 '/')'""'' V - k){n - k\. The SLD 
for this case is A((/)) = t''i''a(f>Y^l=o{9k/ fn~k)\k){k\, 
where gk = 2 [fn-k-i{k + 1)(1 - Sk,n) - fn-kk] 
and fk = ('^!)(sin2 (/.)'= (cos^^)"-''-. The QFI reads 
H{(f)) = tari^ ipj^l^oigl/ fn^k) = 4n. Fock states thus 
enable the optimal unconditional estimation of loss regard- 
less of the actual value of the parameter to be estimated 
(see Fig. [U. This makes an adaptive estimation scheme 
unnecessary (unlike the Gaussian case [1]). Moreover, 
the measurement that has to be performed, obtained by 
projecting onto one-dimensional eigenspaces of A{(j)) H, 
can be implemented only by simple photon-counting. Given 
the recentlyachieved degree of control in this measuring 
technique lllOll and in the high-fidelity engineering of Fock 
states with a small number n < 10 of photons (conditionally 
for running optical fields, and even deterministically in 
microwave cavity or circuit QED) |2;], with n — 2 standing 
as an ideal workpoint, our results might pave the way for an 
experimental verification of the quantum theory of optimal 
estimation and a measurement of the SLD to infer the value 
of such a relevant parameter as the loss factor in dissipative 
channels. While one may argue that in practice it is easier 
to produce and manipulate "classical" fields, i.e. coherent 
states obtained from attenuated laser beams, than nonclassical 
resources such as squeezed (Gaussian) and Fock states, we 
remark that, as shown in Fig. [T] the performance of coherent 
probes for loss estimation is quite far from optimality. In 
particular in the regime of small and intermediate losses 
{(j) < 7r/4), corresponding to routinely available good quality 
fibre channels, the precision achieved by, say, a two-photon 
Fock state would be matched by that of a coherent field 
with much higher mean photon number (e.g. n = 26 for 
(f) = 7r/16, and n = 105 for = 7r/32), an increase in 
energy which may be not worth paying in terms of efficiency 
of the estimation. In fact, in actual implementations it is 
very desirable to have probes of low energy in order not to 
alter the channel significantly 1 1111 and to enable repeatability 
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FIG. 2: (color online) Quantum Fisher information H{(j)) versus the input energy n for different values of the loss parameter (f>. Solid line: 
optimal Gaussian probes fl]; Dotted (magenta) line: qubit-like states \ipl ); Dot-dashed (orange) line: optimal qutrit-like states \il>l ); 
Dot-dot-dashed (cyan) line: optimal quartet-like states |i/;q ). The thin dashed line depicts the ultimate quantum limit, H{(p) = 4n. 



of the input-and-measure scheme. In this respect it appears 
crucial to identify classes of non-Gaussian states which may 
attain the ultimate precision for any value of the input energy, 
especially in the relevant regime of < fi < 1. 

Photonic qubit states. — The simplest candidate probe state 
in the low-energy regime is the superposition of the vacuum 

and the one-photon Fock state, po = l^o )(V'o I' IV'o ) = 
cos 0|O) + e"^ sin 6*1 1), characterized by a mean photon num- 
ber n = sin^ 6. It is rather straightforward to obtain the 
evolved state p^ and to diagonalize it in order to compute 
the SLD. The resulting expression for the QFI is found to 
be independent of the phase tp and given by H^^\(j)) = 
4n[l — (1 — n)cos^(/)]. We notice (see Figs. [T] and |2]) that 
the considered simple example of non-Gaussian superposi- 
tion state (with no free parameter left for optimization) yields 
a significant improvement over the best Gaussian estimation 
for intermediate-high losses, although in the regime of small 
losses and small energies Gaussian states (which in this limit 
are simply squeezed states |llj]) remain better probes. 

Photonic qutrit states. — Next, we consider superpositions 

of the vacuum and the first two Fock States, po — \'4'o /(V'o I' 
with |f/;Q ') = cosQ;|0)+e''' sinasin/3|l)+e*'^ sinacos/3|2). 
Here a can be fixed as a function of (3 and fi, a — 



The evaluation of the SLD involves 



s(2/3)+3^ 

the diagonalization of the 3x3 matrix corresponding to the 
output state. The QFI has to be optimized, for a given n, over 
the phases ji and v and over the weight (3 (the latter ranges 
from /3 = 0, corresponding to a superposition of |0) and |2), 
to /3 = 7r/2, corresponding to the previously considered qubit- 
like state superposition of |0) and 1 1)). Maximization over the 
phases yields ji ^ v ~ n. The optimal (3 can instead be 
found numerically for each fi, 0, and is reported in Fig. |3] 
The resulting optimal QFI is shown in Figs. [T] and |2] The 
photonic qutrit states improve over the qubit-like state and, 
more remarkably, over the optimal Gaussian probes with the 
same mean energy in the whole range of parameters (i.e., for 
any value of the loss). In the limit of vanishing probe energy, 
n — > 0, the optimal Gaussian probe is a purely squeezed vac- 



uum 111] with QFI i/^^) (0) = 4n[l + z2]/[l + 2z(l + n) + z2] 
(where z = tan^ 0), while the optimal qutrit-like state is a 
pure superposition of |0} and |2} (^ = 0) with QFI i/(2)(0) = 



An[l + z'^]/[l + z{2-n + z\ > H^^^ip). In the limit n ^ 1, 
the weight f3 increases up to tt/2 and the qutrit-like state con- 
verges to the optimal Fock state 1 1). 

Relation to de-Gaussifled states. — A very important and 
natural question concerns the nature of such an opti- 
mal qutrit-like state, in particular whether such state can 
be interpreted as a finite-dimensional truncation of an 
infinite-dimensional non-Gaussian state, and how to deter- 
mine the latter. To this aim, we consider de-Gaussified 
photon-subtracted, displaced, and squeezed states pg = 
J\f-^aD{i])S{r)\0){0\S'< {r)D'^ {T])a^ and study their projec- 
tions on the subspace spanned by the vacuum and the first 
two Fock states. The choice for comparison is inspired by 
the fact that the strategies of photon addition and subtraction 
are the current royal avenues to the experimental production 
of optical non-Gaussian resources [12]. Truncation of the 
photon-subtracted displaced squeezed state pp yields the 
state IV'^"'^*''') = co|0) + ci|l) + C2|2), where the coeffi- 
cients Cj = kj/{kl + kl + fc|)^/^, with fco = 77(tanhr + 1), 
ki = k'^ — tanhr, fc2 = 2~^/^fco(fcQ — 3tanhr). The coef- 
ficients kj are functions of the modulus of the displacement rj 
and the real squeezing amplitude r (the relative phase can be 
set to zero in order to maximize the QFI). Remarkably, we find 
that for any -q and r such that < n < 1, the pure states pp 
and pg '^' — |7/)("'5*'"))(7/;("'5*'')| always possess a high- 
fidelity overlap J^ = Tr[p^"^Vo"'^*''^] > 92%. Adding one 
or few further terms in the superpositions quickly rises the fi- 
delity well above 99%. The mean photon number n is a highly 
nonlinear function of -q and r. Thus, in order to visualize 
the maximization of H over these two resource parameters, 
at fixed n, we let -q and r vary in the real space, and depict as a 
shaded region in Fig.[3]the achievable range of (3 (characteriz- 
ing the originally introduced parametrization of fixed-energy 
qutrit states) versus the values of n that are spanned by the 
variation of the parameters. We find that there exist combina- 
tions of values of rj and r for the truncated photon-subtracted 
states such that the angle j3 can take almost any value in the 
range [0,7r/2] for each n S [0, 1]. Superimposing the para- 
metric region with the curves of the optimal /? as a function of 
n yields the maximum QFI among all qutrit states, for differ- 
ent values of the loss parameter 0. Apart from a small range 
of extremely high losses and low energies (in which practi- 







FIG. 3: (color online) Shaded surface: attainable values in the space 
of the coefficients /3 and n of qutrit-like states, as functions of the 
parameters r and rj associated to the truncation of photon- subtracted 
states. Dashed lines: optimal /3 as a function of n, corresponding 
to the maximum quantum Fisher information among all qutrit states, 
for different values of the loss parameter (ranging from 7r/16 to 
7r/2 from bottom to top). 



cally all states, such as non-truncated Gaussian states, qubit 
states, qutrit states, etc.. yield the same optimal performance 
close to the ultimate limit), there always exist values of the 
displacement 77 and of the squeezing r such that the truncated 
photon-subtracted state reproduces exactly the optimal qutrit 
state (pictorially, the dashed lines fall in the attainable shaded 
region in Fig.O. On the other hand, this conclusion does not 
apply to finite truncations of Gaussian states, which are never 
optimal among all pure qutrit states in many ranges of values 
of energies and losses, although they can still perform better 
than the original Gaussian states. This clearly shows that de- 
Gaussification, be it implemented by means of truncation, of 
photon-subtraction, or both, generally enhances the task of es- 
timating loss in bosonic Gaussian channels. It is reasonable to 
conjecture that there exist particular families of non-Gaussian 
states (e.g. non-truncated photon-subtracted states) that repre- 
sent optimal resources for the considered task in the regime of 
low energy, attaining the ultimate quantum limit also for in- 
termediate losses, where superpositions of the first low-lying 
Fock states do not saturate the An scaling of H{<t)). 

Higher-order superpositions. — The previous conclusions 
can be confronted by investigating the effect of adding terms 
of higher order in the superpositions. Consider states of the 

formpo = \'^f^){'^'i\ with IV-o^') = Efl=oC«l?^)> i-e. su- 
perpositions of Fock states up to n = 3. The optimal QFI can 
be obtained by optimizing numerically the complex weights 
c„ for each n, (j) (see Figs. [T] and |2]i. This yields a further im- 
provement over the optimal qubit and qutrit states as well as 
over the Gaussian states. The succession of curves of _ff ('^^ {(/>) 
thus appears to converge to An for fc -^ 00. It is certainly 
true that the best possible performance requires non-Gaussian 
states in every energy range and for any value of the loss. Ob- 
viously, not all non-Gaussian states improve over Gaussian 
ones: For instance, the QFI for cat-like superpositions of co- 
herent states is almost always smaller than that of the optimal 
Gaussian probes, but for n <2 and </> < tt/8. 



Discussion. — Estimating the loss factor of a channel is an 
important theoretical issue of direct practical relevance. In 
the case of purely dissipative bosonic channels, we have 
shown that proper nonclassical non-Gaussian states such as 
Fock states and superpositions thereof need to be employed 
as probes in order to achieve the most precise estimation. 
The optimality of Fock states can be understood on intuitive 
grounds by realizing that the parameter subject to estimation 
is essentially the decay rate of the field energy, and Fock states 
are eigenstates of the energy observable thus having zero en- 
ergy uncertainty in their preparation. 

The present work provides strong support for the need 
of going beyond the Gaussian scenario in applied quantum 
technology and quantum metrology, and motivates further 
research to realize advanced tools of non-Gaussian quantum 
state engineering, manipulation and detection. 

L.A.M. Souza was supported by CNPq-Brasil. 
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